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Abstract — In recent times, we experimentally realized a quite efficient modeling of the shape of 
diffraction-resistant optical beams; thus generating for the first time the so-called Frozen Waves (FW), 
whose longitudinal intensity pattern can be arbitrarily chosen, within a prefixed space interval of the 
propagation axis. Such waves possess a host of potential applications: in medicine, biomedical optics, 
optical tweezers, atom guiding, remote sensing, tractor beams, optical communications or metrology, 
and other topics in photonic areas. In this work, we extend our theory of FWs —which led to beams 
endowed with a static envelope— through a dynamic modeling of the FWs, whose shape is now allowed 
to evolve in time in a predetermined way. And we experimentally create such dynamic FWs in Optics, 
via a computational holographic technique and a spatial light modulator. Experimental results are here 
presented for two cases of dynamic FWs, one of the zeroth and the other of higher order, the last one 
being the most interesting, consisting in a cylindrical surface of light whose geometry changes in space 
and time. 

PACS nos.: 42.25.Bs; 42.25.Fx; 41.20.Jb; 46.40.Cd; 41.85.-p; 46.40.Cd; 42.20.Ht; 42.40.Jv, 42.30Lr. 

1. Introduction — The “Frozen Waves" (FWs)[Tl[^|3] are very peculiar solutions to the wave equations 
consisting in non-diffracting beams whose longitudinal intensity pattern can be freely chosen within a 
prefixed space interval 0 < z < L oi the propagation axis. For such beams even a certain degree of control 
on their transverse intensity pattern is possible. This endows the FW-type beams with a great deal of 
potential applications, since the concrete possibility of modeling their spatial intensity shape adds to the 
exceptional characteristic of all the “Localized Waves", or Non-Diffracting Waves (NDW), of resisting the 
effects of diffraction in |5J E] . 

The theoretical development of FWs can be found in Refs. nnamiii, besides in the above-quoted 
reviews or books: We call it, briefly, the FW method. Synthesizing, a FW is a superposition of co¬ 
propagating and equal-frequency Bessel beams (of zeroth or higher order) 

N 

n=-N 

with = ^/[Jv{-)\max, where [Jv{-)]max is the maximum value of the Bessel function of the first 
kind Jv{-)', kpn and kzn are transverse and longitudinal wave numbers of the n-th Bessel beam in the 
superposition, satisfying the known relation —The FW method requires the longitudinal 

wavenumbers to be chosen so that kzn = Q + 2'Kn/L, quantity Q being a constant to be selected, in the 
case for instance of zero order FWs, according to the desired spot-size Apo, besides for constructing solely 
(not evanescent) beams propagating in the positive direction. (Let us recall that a partial control on the 
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beam transverse intensity too can be obtained by a suitable choice of Q, and also by having recourse to 
higher order Bessel beams). 

An important point is that the FW method [TJ |3J |7] shows that, once it has been arbitrarily chosen 
the desired intensity longitudinal pattern \F{z)\'^ in the spatial interval 0 < z < L, then the coefficients 


in superposition y have to be A„ = (1/L) F"(z)e“*T?”^dz. 

This is enough to obtain Diffraction-Resistant beams with the required longitudinal intensity pattern, 
concentrated: (i) either along the propagation axis (p = 0), when zero-order [y = 0) Bessel beams are 
chosen in solution 0 ; it is then possible to determine the spot radius of the resulting beam from the 
parameter Q via the relation Q = [yp' jc? — 2.4^/Apq)^/^; (ii) or over a cylindrical surface, if a positive 
integer is adopted for z/ in Q ; in which case, the cylinder radius can be approximately evaluated through 
the equation \{AIy{py/jP — Q^)] = 0. 


P=Po 

It should be stressed that optical FWs have been concretely generated via computational holography: 
see Refs. BE]; without forgetting that ultrasonic FWs have been reproduced also in Acoustics, at least 
by simulated experiments ITU]. The experimental production of the FWs opened the door to a great deal 
of applications in a number of different areas: like medicine, biomedical optics, optical tweezers, atom 
guiding, remote sensing, tractor beams, optical communications or metrology, and so on. 

Before going on, let us recall that FWs belong to the set of the subluminal NDWs, corresponding 
to the very particular case of zero Deak-velocitv fTTl I12| : They possess, indeed, a static envelope (within 
which only the carrier wave propagates) [TJ B [3|- appears convenient, however, investigating, and 
constructing dynamic FWs, for all applications in which a traveling FW appears to be desirable and 
useful. Aim of this work is just extending the FW method towards a dynamic modeling of the non¬ 
diffracting beams, whose shape can slightly evolve in time in a predetermined way. That is, now we 
want our resulting non-diffracting beams to be endowed with a spatial shape which changes in time, as 
requested by the particular application considered; and to produce experimentally such dynamic FWs in 
Optics, via a computational holographic techniquejHJ [9l 1181 114| in which computer generated holograms 
of the FWs are created, as specified below, and then reproduced by use of a spatial light modulator 
(SLM). [Such a resulting wave, that corresponds to an intensity pattern no longer static, and slightly 
loses therefore its monochromaticity, could not be called a beam, strictly speaking...: However, if we 
consider the time evolution of its shape much slower than the period associated to the relevant frequency, 
our approach ends constituting an excellent approximation]. 

In equation Q we need inserting the time dependence, first, of the coefficients A„, and, second, of 
the parameter Q which enters the definition of the wavenumbers. Once it has been chosen the desired 
longitudinal intensity |F’(z,t)p, one can exploit the way the said temporal variations are implemented, 
so as to get the wished dynamical control on the longitudinal shape of the FW (and at a certain extent 
on its transverse shape too). 

2. Theory of the “Dynamic Frozen Waves" — Let us consider the time-dependent solution 


N 

An{t)Jj, {kpn{t) p) exp{iiy(j))exp{ikzn{t)z) exp{-iujt) (2) 

-N 


where one may recall that any expression of the type Ju{kpnp) exp(ifc 2 „z) exp(—iwt), with kp„ + k‘^„ = 
uPje?, is a solution to the wave equation (V^ — XlePdfy^ = 0. As already mentioned, in the case of 
slow temporal variations, we are able to obtain in the dynamical regime the longitudinal pattern A(z, t) 
and either the spot radius Apo{t) in the case of a zeroth-order FW; or the radius poit) of the intensity 
cylindrical surface in the case of a higher order FW. 

Namely: To obtain a dynamic FW, instead of F{z) we have to consider for 0 < z < L a time- 
dependent function F{z,t). Then, considering superposition Q, one can determine the coefficients A„(t): 


2 
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An{t) 


- F(z,t)e-*T-dz, 


(3) 


with the time-variation of the longitudinal or transverse wave numbers = a/ (2tt/X)‘^ — and 

kzn(t) = Q{t) + 27rn/L depending on Q(t)-, which in its turn is linked either to the spot size ApQ{t) 


Q{t) = VP-(2.4/Apo(t))" 


in case of zero order FWs; or to the cylinder radius po (t) 


( 4 ) 


[{d/dp)J^{py^LOycP 


QHm 


P=Po(t) 


= 0 


( 5 ) 


in case of higher order FWs. 

We anticipated that a slow variation of F{z,t) and Apo(i) (or po{t)), guarantees Eq.([^ to be an 
approximate solution to the wave equation. This can be shown by noting that for each term of the 
superposition Q we have 


-^[An{t)J^{kpn{t)p)exp{ikznit)z)exp{-iwt)] « An{t)J^{kpn{t)p) exp{ikznit)z)-^ exp{-iu}t). (6) 

The above approximation is valid due to the said slow variation, which implies that 

<< \ujA„{t)Ji,{kpn{t)p)\. All this guarantees that Eq.(|^ is a good ap¬ 
proximate solution to the wave equation. 

3. Generation of Dynamic Frozen Waves (DFW).' Experiments and results — As already mentioned, 
the experimental holographic method for the generation of static FWs can be found in Refs. mi. Our 
present setup is shown in Fig[^ It is composed by an Argon ion laser (of 514,5 nm), which is expanded 
and collimated (Exp) in a nematic SLM device (LC R1080 SLM, HoloEye Photonics AG). We perform 
the amplitude modulation of a computer generated hologram (CGH) on the SLM by polarizer Pol and 
analyzer Anl, with angles 0° and 90°, respectively, measured with regard to the input SLM axis. The 4-f 
spatial filtering svstemfTHIHlin] is used for the FW generation, the SLM (GGH) being placed at the input 
plane (focus of the first lens); while a spatial filtering mask (SF, band-pass circular pupil), located at the 
Fourier plane, selects and transmits the light from the initial spectrum, generating the encoded signal at 
the output plane of the setup. As a result, we get a propagation of the FW, whose intensity is registered 
with a GGD camera that can be displaced all along the distance 0 < z < L. 

Specifically, in the case of DFWs we discretized the relevant time interval (of 10 s, for instance) into m 
parts; and we obtained, for each value Tm of time, the coefficients of the series and consequently the 
complex function 4t(p, z,t = Tm) describing the FW at that instant of time. On knowing the field value 
at each time instant, m computer generated holograms (GGHs) are created. Then, by using a Spatial 
Light Modulator (SLM) linked to a microcomputer, the holograms are reproduced at a specific rate on 
the modulator display to get a reconstruction of the DFW: That is, of position and longitudinal (and 
transverse) shape, varying with time, of the dynamic frozen wave (see FigQ. The LG-R1080 modulator 
possesses a rate of 70 Hz, enough for producing 15 frames per second. 

The z-axis too is divided into j parts Zj. The data acquisition is then performed leaving the GGD 
camera in the fixed position zq = 0 (the beginning of the propagation axis), where all frames Tm repro¬ 
duced by the SLM are captured; the process is repeated in the subsequent positions till the j-th one, as 
in Fig|^ 
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Figure 1: (color online) The experimental setup used for generating a dynamic Frozen Wave is essentially 
composed by an Argon ion laser beam, FLF being a beam collimator, Pol a polarizer, BS a beam splitter, 
SLM a nematic Spatial Light Modulator, and Anl an analyzer; while the L’s are lenses, SF is a spatial 
filtering mask, and CCD is a camera for images/data acquisition. 


3.1. First example (for a zeroth-order DFW): - Let us consider, in the space interval 0 < 0 < L, an 
exponentially rising barrier of the type 


F(z,t) 


exp(Arzt), for Zi < z < Zf 
0, elsewhere 


( 7 ) 


with Zi = zq + vt and Zf = AD + zq — vt, quantity AD = Zf — Zi being the width of the considered 
DFW at the initial time instant, and K~^ = 1 cm s. The coefficients of our series, at each instant t of 
time, are given by Eq.([^ when inserting into it the function ([^. Afterwards, the FW behavior in the 
considered time interval (e.g.. At = 10 s) is obtained by putting such values A„(<) into solution ([^. 
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Figure 2: Scheme of the experimental data acquisition. The z-axis is divided into j parts the data 
acquisition is performed leaving the CCD camera in the fixed position zq = 0, where all frames 
reproduced by the SLM are captured. The process is repeated in the subsequent positions till the j-th 
one. 


Let us adopt the values: Q = 0.99997 w/c, and L = 40 cm, with = 20 superpositions; a step 
width AD = 4 cm; initial position zq = 0.18 cm, and finally a constant v = 1.3 cm/s. The initial 
and final positions move with velocity —v and v, respectively. We then get a step profile which, as time 
elapses, modifies into the longitudinal exponential ramp expressed by Eq.Q. More specifically, let us 


4 




















divide the time interval At = 10 s into m = 150 parts: From Eqs.(|^ and ([^ we obtain z,t = Tm), 
and for each time value an amplitude hologram is generated, with a total of 150 holograms (frames). 

In the present case, the z-axis was divided into 81 parts, and 143 frames were captured at each 
position Z{j). This implies an error in the temporal reconstruction of the holograms (error due to 
the intrinsic limitations of the microcomputer hardware and of the software itself: think of the time 
spent by Matlab between command and execution) that can be evaluated to be of about 7%. In Figsj^ 
the longitudinal intensity profile for different time instants is shown, together with its corresponding 
orthogonal projections (while online, in “Media-1" the 3-D profile is shown in a dynamic way). 
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Figure 3: (color online) Time evolution (for At = 10 s) of the longitudinal intensity pattern for the FW 
corresponing to the longitudinal intensity profile given by an exponentially rising ramp like ([^. In the 
insets, the corresponding orthogonal projection are shown. See also Media-1 online. 

3.2. Second example (for a first-order DFW): - We have presented above our experimental results 
for a dynamic FW of the zeroth order. The most intriguing cases are, however, the ones dealing 
with higher-order dynamic FWs. For the case of higher-order (static) FWs, in previous papers we 
theoretically constructed electromagnetic or acoustic static envelopes with the shape of donuts, or of 
cylindrical surfaces [m I213 [H] . Going on to higher-order DFWs, striking results appear: We present in 
our last Figure the experimental results referring to a first-order DFW, consisting in a cylindrical surface 
of light whose geometry changes in space and in time. 

Our theoretical choice for a first-order DFW is as follows: We may theoretically construct a higher- 
order DFW, characterized by a duration T, whose behavior may be divided into three steps, with a 
duration r/3 each, so that Ti = T/3 and T 2 = 2T/3; and with a deformation speed (in each step) given 
by the alteration length divided by T/3. Namely, let us devise a FW initially consisting in a short 
cylindrical surface of length Ali and radius pi, such that: (i) in the first step, the cylindrical surface is 
required to stretch, keeping its radius, till reaching the length Alf, (ii) in the second step, the cylindrical 
surface keeps its length, but is requested to reduce continuously its radius, till the value pf\ (iii) in the 
third step, the cylindrical surface —while keeping its radius pf — is required to split into two cylinders, 
which become shorter till reaching the length Alf each. 


5 















Let us choose (for all steps) the simple function F{z, t) defined as zero everywhere (within 0 < z < L) 
except 

{ 1 for zi < z < Z 2 

( 8 ) 

1 for Z 3 < Z < Z 4 , 

depending on time through the values of zi, Z 2 , Z 3 and Z 4 . Its radius will change in time according to 

Q = Q{t). 

To be more specific: 

(i) first movement, 0 < t < Ti = T /3: During this phase, the cylinder radius po = Pi does not change, 
and the corresponding value of Q remains fixed, according Eq. The z values are Zi = Zii — Uit; 
Z 2 = Z 3 = zii +AZi /2 = z; Z 4 = (zii +A/i)+nit, that is, zi and Z 4 change in time, while Z 2 and Z 3 remain 
fixed. The initial values of zi and Z 4 , that we call zii and Z 4 i = zi^ + /Sli, are still to be fixed; together 
with the final value zi/ of zi, with zi(Ti) = zi/ < zii, which also determines the final value of Z 4 , so that 
Z 4 / = Z 4 (ri) = (zii + ISdi) + (zii — zif). The final length of the cylinder will he Alf = Z 4 / — zif, while 
the speed vi will be vi = (zu — zi/)/Ti = 3(zii — zi/)/T. 

(ii) second movement, Ti < t < T 2 = 2T /3: The z values are the same assumed at time f = Ti in the 
previous step, that is, zi = zi/; Z 2 = Z 3 = z; Z 4 = Z 4 /; but quantity Q will vary according to Eq.(|^, 
where po(t) = pi — V 2 (t — Ti). The final value of the radius, ^ 0 (^ 2 ) = P/j is still to be fixed, our request 
being merely that pf < pi. The speed V 2 will be V 2 = {pi — pf)/{T 2 — Ti) = 3{pi — pf)/T. 

(iii) third movement, T 2 < t < T 3 = T: Quantity Q keeps constant and equal to the final value it 
acquired in the previous step. The z values will however change as follows: zi = zi/; Z 2 = z —n 3 (t —T 2 ); 
Z 3 = z + n 3 (t — T 2 ); Z 4 = Z 4 f. The final value of Z 2 (and of Z 3 ) remains to be fixed, for instance through 
a choice for the final length (AZy) of the two smaller cylinders: In fact, it will be Z 2 (T') = zi/ + AZ/; and, 
for the final value of Z 3 , one will have Z 3 (T) = Z 4 / — Alf. The speed is ^3 = (z — Z 2 {T))/{T 3 — T 2 ) = 
3(z-Z2(T))/T. 

This specifies all what is needed for the experimental creation of the devised higher-order (first order) 
DEW. 

Its experimental production: — The experimental procedure is the one illustrated in the previous 
sections 2 and 3. For the definition interval of the DEW, that is, 0 < z < L, it has been chosen the value 
L = 40 cm, being the other parameters chosen as: T = 10s, zii = 0.18 cm, zif = 0.05 cm, Ak = 0.004 
cm. Pi = 35 ^m and pf = 25 /rm. With this we have vi = 0.039m/s, v2 = 3/im/s and v3 = 0.03m/s. As 
to Eq.Q, the value fV = 12 has been adopted. 

In Figs|^the experimental results are shown of the time evolution, during At = 10 s, of the inten¬ 
sity pattern for our first-order dynamic FW; as well as the orthogonal projections of the corresponding 
cylindrical surfaces. The final radius of the cylindrical surface has been chosen to be p/ = 25 pm, and 
its initial radius results to be pi = 35 pm. Noticeably, as time elapses, the cylindrical surface stretches 
its length (while narrowing its radius), and afterwards splits into two short cylinders: see also Media-2 
online. 

4-. Conclusions — The Non-Diffracting Waves (NDW) are solutions to the (linear) wave equation 
which travel well confined (localized), in a single direction, and with other interesting properties like 
self-healing. The Frozen Waves (FW) are NDWs with zero peak-velocity, i.e. with a static envelope, 
whose longitudinal intensity pattern can be freely modeled —both theoretically and experimentally— 
within a prefixed space interval of the propagation axis. In this paper we present, first, the theory of 
“dynamic FWs" (DFW) whose shape is now allowed to slightly evolve in time in a predetermined way, and, 
second, their experimental generation by holographic techniques (namely, by a holographic setup along 
with an optical reconstruction of dynamic holograms CGHs by using a LC-SLM). The most interesting 
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Figure 4: (color online) Time evolution, during At = 10 s, of the intensity pattern for our first-order 
dynamic FW; and orthogonal projections of the corresponding cylindrical surfaces: experimental results. 
The final radius of the cylindrical surface has been chosen to be p/ = 25 p, and its initial radius results 
to be pi = 35 p. Noticeably, as time elapses, the cylindrical surface stretches its length (while narrowing 
its radius), and afterwards splits into two shorter cylinders: See also Media-2 online. 


experimental results, depicted in our last Figure (and in movies online), refer to the higher-order FW 
consisting in a cylindrical surface of light whose radius changes in space and time. We show that, as 
time elapses, the cylindrical surface stretches its length (while narrowing its radius), and afterwards splits 
into two shorter cylinders. The experimental results are in agreement with the theoretical predictions, 
and open exciting possibilities of generating many further potentially interesting DFWs for scientific and 
technological applications. 
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